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Abstract. We consider the long time behavior of solutions to the magneto- 
hydrodynamics equations in two and three spatial dimensions. It is shown 
that in the absence of magnetic diffusion, if strong bounded solutions were to 
exist their energy cannot present any asymptotic oscillatory behavior, the dif- 
fusivity of the velocity is enough to prevent such oscillations. When magnetic 
diffusion is present and the data is only in L 2 , it is shown that the solutions 
decay to zero without a rate, and this nonuniform decay is optimal. 



1. Introduction 

We consider the Magnetohydrodynamics equations (MHD) in two and three 
dimensions. We deal with questions regarding long time behavior of solutions to 
the MHD with and without magnetic diffusion. The MHD equations model the 
interactions between a magnetic field and a viscous incompressible fluid of moving 
electrically charged particles. 

In non-dimensional form the equations can be expressed by 

+ (u ■ V)u - S(B ■ V)B + V(P + | \B\ 2 ) = J-Au, 
at 2 Re 

+ (u ■ V)B - (B ■ V)u = SAB, 

V • u = 0, V • B = 0, 

u(x,0) = uq(x), B(x,0) = B (x), 



where u = u(x, t) = t), u 2 (x, t), . . . , u n (x, t)), B = B(x,t) and P = P(x, t) 

denote the unknown velocity, the magnetic field and pressure of the fluid at a point 

\B\ 2 

(x,t) £ K™ x R + , respectively. The term J-J- denotes magnetic pressure. The 
positive constants appearing in the equations are Re, the Reynolds number, Rm, 
the magnetic Reynolds number, and S = M 2 /(ReRm), where M is the Hartman 
number. For the sake of notational simplicity, and with minor loss of generality, we 
set all these constants equal to one. After rescaling u and B, let p = P + \B\ 2 
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denote the total pressure, Equation (1.1) can be rewritten as 
d 



at 



u+(u- V)u - (B ■ V)-B + Vp = Aw, 

rj 

— B + (u ■ V)B - (B • V)u = <5AB, 



V • u = 0, V • B = 0, 
u(x,0) = Uo{x), B(x,0) = Bq(x). 



(1.2) 

(1.3) 
(1.4) 
(1.5) 



The initial data (uq(x), Bq(x)) will be chosen below in appropriate spaces. Deriva- 



tions of these equations can be found in [Cha81, Cow76, LLP84|. 

Many authors have studied MHD equations from the point of view of existence 
and long time behavior. Without making a complete list of all authors we would like 
to mention some of the relevant literature. In particular important results on exis- 
tence were obtained, among others, in HKoz87| |HX05bj |HX054 |Wu02j . In the direc- 
tion of decay interesting results can be found in the papers pS89j [Kim02|, |SSS96| . 
The methods used for decay in [MS8E] were based on Fourier Splitting [ Sch85| . The 
paper [Kim02| uses Gevrey regularity and it is based on ideas developed by Foias 
and Temam. Similar methods were used for the Navicr-Stokes equations and can 
be found in [|OT00 



Several questions will be addressed regarding the long time behavior of the solu- 
tions. In this regard, in the absence of magnetic diffusion, we are going to analyze 
the following problem. 

Is the diffusion introduced by the velocity alone sufficient to prevent compensatory 
oscillations? Specifically, simple calculations shows that the following energy in- 
equality holds when (5 = 0, 

±(\\ u (t)\\l + \\B(t)\\l) <-2\\Vu(s)\\l 
This shows that the combined energy decays, but allows the possibility of separate 



oscillations in 



and || .B || 2 that could compensate each other. In this paper we 



show that in the case that there were to exist strong solutions which are bounded 
such oscillations can never occur. Specifically it is shown that 

Theorem 1. Let n = 3. Suppose that there exists (u,B) a strong bounded solu- 
tion to the MHD equations with 5 — 0, magnetic field bounded in L°° , and data 
{uq, Bo) e (L 1 fl L 2 x i 2 fl L°°). Then 



\\u(t)\\ 



0. 



\B(t)\\ 



M, 



as t — > oo, where M is some positive constant. 



The previous theorem shows that the diffusion in the velocity is sufficient to 
prevent compensate oscillations between the two energies. 

The next question we addressed is in regards to decay of solutions to MHD 
equations with diffusion both in the velocity and the magnetic field. Here it is 
supposed that the data is only in L 2 space. In this case it is shown that the energies 
of the velocity and the magnetic field decay to zero without a rate. Moreover, it is 
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shown that this is optimal, that is, that cannot be a uniform rate for the energy of 
the solutions with data exclusively in L 2 . We show that 

Theorem 2. Let n = 2, 3. If (it, B) is a weak solution to the MHD equations with 
S > and data (u ,B ) € (L 2 (R n )) , then 



Urn (\\u(t)\\ 2 2 + \\B(t)\\l) =0. 



The proof of this theorem is given first formally. To make the result rigorous we 
apply the method to smooth approximations and then pass to the limit. 
With regard to the optimality of this last result we show that 

Theorem 3. There exist no functions G(t,f3) and H(t,~/) with the following two 
properties. If (u,B) is a solution to equations ( |-?.l[ j-( [77^ j with S > and data 
(u ,B ) G (L 2 (]R n )) 2 , n = 2,3, then 

i) ||u(t)|| 2 <G(t,||«o|| 2 ), 

ii) lim t _> G(t,0) = 0, for all /3 > 0. 

The last part of the paper focuses on extending Kato's pioneering work on L p 



decay for Navicr-Stokcs equations [Kat84] to the MHD equations with magnetic 



diffusion. In particular we note that a simple modification of Kato's work |Kat84| 
yields equivalent results for the MHD equations, then combined with our result on 
non-uniform decay in L 2 gives a slight improvement of the decay rates. 

Corollary 4. There is X > such that for ||uo|| 2 < A the global solution of the 



equation (6.1) for q > to, and for 2 < r < q 

hm t^r ||(u,fl)|| r = 



2. Notation 



We denote 



Cq° ct = Cq^IW 1 ) : space of infinitely diffcrcntiable functions with compact 
support and divergence free, 
L 2 = L 2 (R n ) : completion of C™ a in the L 2 -nonn ||-|| 2 , 
Hp = _ff*(]R n ) : completion of C^ a in the homogeneous 7? 1 -norm ||V-|| 2 . 

The Fourier transform of cp will be denoted by ^{f} = <p = e~ l ^ x ip(x) dx and 
its inverse transform by .^^ 1 {ip) = <p = / e lx ^tp(£)d£. Also, 




L p (a,b;L") = if : (a, b) x ^ : \\f\\ LP{aMLq) = / ||/(r)||J dr < oo 



The notation || -\\ LP , q will be used for the norm of L p (0, oo; L q ), and (/, (?) = J fg dx 
for the inner product in L 2 . Various constants are simply denoted by C. 
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3. Preliminary Results 



In this section auxiliary results which will be needed in the sequel are obtained. 
We consider the MHD equations with no magnetic diffusion (S = 0). Some of the 
results presented are standard and their proofs arc included in the appendix for 
completeness. The results below follow ideas of JORS97 ]. 

We start with an estimate for the Fourier transform of the velocity. 



Proposition 1. Let (u,B) be a mild solution to the MHD equations (l.i)-(1.5) 
Assume the initial data i*o, Bq is in L 1 (M 3 ) n L 2 (R 3 

1 



The 



\u(t)\ <C 1 + 



where C is a constant. 



Proof. See Appendix, Proposition 0. □ 
The proofs in this section are formal. To make them rigorous is suffices to apply 
them to approximations using retarded mollifications such as the ones constructed 
by ]CKN82| , |Ler34| , and pX05b{ for the MHD equations and then pass to the 
limit. For details regarding our proofs see ]ORS97 were the procedure has been 
done for the solutions to the Navier-Stokes equations. The extension to MHD is 
straightforward. 

We first recall a standard energy inequality 



Ht)\\l + \\B(t)\\ 



| Vm(s) || 2 ds = 1 1 XtQ 1 1 2 



|-B || 2 



This follows easily by multiplying fll.2| ) by u, ( jl.3| ) by B, adding the equations, and 
then integrating in space and time. 

The next proposition gives a generalized energy inequality 



Proposition 2. Assume (u,B) is a weak solution of Equations (l.i)-(l.L). For 
E{t) e C^R; R + ) with E(t) > and tp € C^R; C^OL 2 ), the weak solution satisfies 



E{t)U*u{t)\\l=E{s)U*u{s)\y 2 + / E'(t)U*u{t)\\ z 2 dr 

J S 

+ 2^ Eir)^' *u(t),iP*u{t)) -\\Vtp*u(T)\\fj dr 
2 f E(t)((u-'Vu(t),'iP*iIj*u(t)) 

+ (B ■ V-B(r), ip * ij) * u{t))) dr 



(3.1) 



Proof. Multiply Equation (1.2) by E{t)ip * ip * u (t) an d integrate by parts to get 
f t (E(t) ||V * u(t)\\l) = E'(t) U * u(t)\\ 2 + 2E(t)W * * «(*)> - IIVV * u(t)f 2 } 

- 2E{t){(u ■ Vu,ip * ip * u(t)) - (B ■ V.B, ip*ij)* u(t))} 

Integrating the preceding equation in the interval (s,t) yields Equation ( |3.1| ). □ 
Corollaries || and || follow as an easy consequence. 
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Corollary 5. Let (u,B) be a weak solution of (l.i )- §U^ . Let tp e L 2 (]R 3 ), then 
||0*«(t)||a< e A{t - s} <p*u{s) \2 f (\(u ■ Vu,e 2A{t ~ T) £ 2 * u{t))\ 

^ J S 

+ \(B-VB, e 2A(t ~ T V 2 * u(r)) I ) dr. (3.2) 



Proof. Choose E(t) = 1 and ^p(r) as 

^(r) = ^" 1 {e- | « |2(t+I '- T V(-e)}, V > 

in Eq. (^TTJ) . Then ^(r) * w(r) = e^+'^V * «(r) = / e At - t+, '~ T+s ^(x)u(6) d5, 
and 

(ip'ir) * w(r), -0(r) * it(r)) — ||V?/>(t) * u(r)|| 2 

= -(A(e A(t+T '- r) (^) * u( T ), e A (*+"- r) (^ * u(r)) - V(e A (*+"- T) ^) * u ( r ) 
= 0. 



Hence we have from (3.1) 
||e A "(^*u(t)|| 2 < 



°'tp * u(s) 



+ 2 y Vii,e 2A < t+ ^ T ^ 2 «(T)) 

+ (B- VB,e 2A(t+ '- T ^ 2 *ii(r)) 



r/r 



Let fy — > in the preceding equation to obtain (3.2). □ 
Corollary 6. Let E{t) G C 1 (IL + ;K) and ip G C 1 (0, oo; L°° n L 2 ). T/ien a weak 



solution of Equations (LI )-( j-?.q ) satisfies 

E(t)U(t)u(t)\\l <E(s)U(s)u(s)\\l+ I E\t)U(t)u{t)\\ 2 2 dr 



+ 2 / B(r) (^(r)«(r),V(T)ti(r))- f^(r)u(r) , dr 



(3.3) 



£(t)((^{u ■ Vu{t)},^ 2 u(t)) 

+ (^{B -VB(t)}J' 2 u(t))^ dr, 
for almost all s > cmd all t > s. 
Proof. Apply Planchcrel's theorem to (|3.l[). 



□ 



4. Non-uniform decay of solutions to the MHD equations with no 

MAGNETIC DIFFUSION 

In this section it is shown that if there were to exist strong bounded solutions in 
3D, they can not have compensatory oscillations. We analyze separately the energy 
of the high and low frequencies of the solutions. The main tool for the analysis of 
the high frequency is Fourier Splitting, see [Sch85|. 
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We first establish Theorem [l] of the introduction, which we recall for complete- 



Theorem 1. Let n = 3. Suppose that there exists (u,B) a strong bounded solu- 
tion to the MHD equations with 6 = 0, magnetic field bounded in L°° , and data 
(uo, B ) G (L 1 n L 2 x L 2 n L°°). Then 

|Kt)|| a -0, and \\B(t)\\ 2 ^M, 

as t — > oo, where M is some positive constant. 

Proof. Split the velocity of the solution into low and high frequency parts 
IK*)ll2 = l|w(i)|| 2 < \\vu\\ 2 + \\(l-v)ii\\ 2 , 

where ip is a function in Fourier space to be chosen appropriately, to emphasize the 
low and high frequency of u. 

Low frequency Decay. Set ip(£) = e - '^ using the result of Corollary || and 
Plancherel theorem, 



< 



e-^^cpu 



(s) 2 + 2 J s («'Vti,e 2AM ^ 2 *«} 



[B- VB,e 2A ^- T ^ 2 *u) 



< 



-\Z\ 2 (t-s) 



ipu(s) 



2 f (ip 2 * B ■ V-B, e 2A ' t_r 

J S 



(ip 2 *u- Vti,e 2A(t - T) ii) 
,1t 



dr 



Clearly the first term on the right hand side satisfies 



-I?| 2 (*-s) 

To bound the third term note first that 



limsup||e 141 ^ s> ipu(s)\ 

t — >oo 



0. 



(4.1) 



\<P * 



B - VB \\l = H J \$ 2 *{B- V)B ] \ 2 dx < CJ2 (j2\\ d ^ * ^^h) 

j 3 V i I 

<C||B!|^^||^ 2 || 2 ^ , 



where C is a positive constant. 
Thus 



{(p 2 * B ■ VB,e 2A ^u)\ < C \\B(r)\\l \\u(r)\\ 2 \\di 



Since ip 2 = (47rf)- 3 e - |j:|2 / 2t , it follows that 

\\d it p 2 f % <C f (t- A \x\e-^^) 2 dxKCr 11 / 2 



R 3 



-i*r/2t 



' -dx^cr 11 ' 2 . 



(4.2) 



NON-UNIFORM DECAY OF MHD EQUATIONS 



7 



Hence, 



(cp 2 * B ■ VB, e 2A (*- r M < C \\B(t)\\ 2 \\u(t)\\ 2 r 11 ' 2 < Ct" 11 / 2 



since the L 2 -norm of u and B are bounded by the initial data. 



Similarly, | (p 2 * u ■ Vu, e^'^u) \ < Ct' 11 / 2 . Hence by (JO]), Q, and the 
last two inequalities 

||<^u(i)|| 2 — > as t — > oo. 
High frequency Decay. We will show that 

lim ||(l-^)u|| 2 <e (4.3) 

for all £o > 0. The Fourier splitting method will be used. Let y(e) = {£ : l£l < 
G(e)}, a neighborhood of the origin, were G will be specified below. Set ip = 1 — p, 
where tp is given above. Then ijj' = |£| (p, and Corollary (^|) yields 

E(t)\\(l-<p)u(t)\\l<E( S )\\(l-p)u( S )\\ 2 2 + f E'{t) I |(1 -p)u{r)\ 2 didr 



+ / E'(r) 



|(1 -p)u(r) | 2 d£dr 



2 / B(r)||e(l- V )fl(r)||* dr 



+ 2 / £(t)(|£|V(tHt),(1-^(t)Ht)) dr 
+ 2 / £(t) (iTVB, (1 - <^) 2 w(t)) c 



+ 2 / E(t) (u-Vu,(l~ip) 2 u(T)} 



(4.4) 

The terms in the second and third row are bounded by 

(E'(r) - 2E{r)G 2 (e)) [ |(1 - p)u(r)\ 2 d^dr. 

Choose E(t) = e £t and G(e) = hence E'(t) - 2E(t)G 2 (e) = 0, thus the above 

integral vanishes. 



Divide Equation (4.4) by E(t), 
||(1 -<p)u(t)\\l< ||y 11(1 -<P)u{s)\\l + 



E(t) 



E'{t) \(l-p)u(r)\ z d£dr 

hie) 



E{t) 
2 

W) 



E(tM\ 2 V (t)u(t),(1- V (t))u{t)) dr 



E{t) (wVw,(1-^) 2 «(t)) 



dr 



E(t) 



I{t) + II(t) + III(t) + IV{t) + V(t). 



E(t) (B^VB,(1~p) 2 u(t)) 
(4.5) 



dr 
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Estimate for I(t): 

I(t) = §A 11(1 - <PM*)\\1 < z e[s ~ t] < (hoWl+WBoWfje^-* < Ce«'-*l 



Hence 



lim I(t) = lim §4 IK 1 " PM^Wl = 0- 



t"oc^ y "' t^oo E(t) " V " T-/ "V"/ H2 

Estimate for II (t): By Proposition |l| 

Jx(s) Jx(e)\ \t,\J Jx(e)\ |£| J 

Since (1 - ipf < 1, 

H (t) = -J^f ' E\t) I \{l- V )u{r)\ 2 d^dr 



E{t) 

f I \Hr)\ 2 d£dT<CeV 2 = e . 

E(t) Js J x {e) 

Estimate for III(t): Observe that < ip - ip 2 < 1, and E(t) < E(t) for r < t, 
hence 



IH{t) = E(r) (|£| 2 <p{t)u(t), (1 - ¥>(r))fi(r)) dr 

E(t) I \i\ 2 {p~V 2 )\u{T)\ 2 d^dr<C [ \\Vu(r)\\ldr, 

Jm, 3 Js 

Since / °° || Vw|| 2 < 00 it follows that 

lim III{t) < C lim lim / ||V«(r)||o dr = 

t — >00 S — >00 t — >00 J 

Estimate for IV(t): Set ( = - (1 - ^) 2 }- This function is essentially 

the heat kernel. Note that (u ■ Vu, u) = 0, hence IV (t) can be estimated as follows 



E(t) 

< 



E{t) 



2 '* 



IV ^ = EmJ E(t)\(u-Vu,(*u}\ dr<C J ||CII 6/5 II«-V U || 3/2 || U || 2 
< C j ||C|| 6/5 ||u|| fl ||Vu|| 2 dr < C f ||C|| 6/5 \\Vu\\l dr 

J s J s 

<Cp^*<££\\Vu\\l*r. 

Thus 

lim IV(t) < lim lim -^77 / ll Vu ll2 dr = °- 

t—>00 S— *OC t-^OO S ' J g 

Estimate for V{t): 



d.T 
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V(t) 



< 



< 



< 



< 



E(t) 
C 

m 

c 

E(t) 

C 
E(t) 

C 
E(t) 



E(t 



E(t 
E(t 
f E(r 

J s 

E(t 



Here is the only place 
we use that B £ L°° 



(B ■ VB, (1 - ip) 2 u(t)) 



dr < 



C 

W) 



/V)(e/|^ 



d£ dr 



\B-B\ Z d£ 



r \ 1/2 

J \B(r)f dej HV«(r)|| 2 dr 
BW^WB^WjVuir^dr 



|Vu(r)r d£ dr 



V, 



||Vii(T)|| 2 dr, 

where we need that the magnetic field is bounded. Specifically 
R 3 x Since B G L 2 we get the above bound. Recall that 



E(t) 



- st 



v(t)< 



hence 

C 
E(r) 

C 



J*E(t)\\Vu(T)\\ 2 dT<^(^p) ' ^\\Vu(r)\\ 2 2 dr 



1/2 



|Vix(t)|| 2 dr 



1/2 



which as before tends to zero as t and s goes to infinity. 

Combining the estimates I(t) — V(t) yields In summary, we have showed that 



lim ||(1 - <p)u(t)\\ 2 < £ . 

£—►00 



(4.6) 



Since £o is arbitrary and positive, combining (4.1) and (4.6) yields 

lim Kt)|| 2 =0. 
To obtain the limit of ||^|| 2 proceed as follows. Set 

<t>(t) = \Ht)\\t+\m)\\l- 

Given that <fi(t) > and is decreasing, there exists a constant M such that 4>(t) 
M ast-^r oo. Since ||u(i)|| 2 — > 0, it follows that 



\B(t)\\ 



M as i-> oo. 



This completes the proof. 



□ 



5. MHD EQUATIONS WITH DIFFUSION 

In this section it is shown that if the data is only in L 2 then the solution decays 
without a rate. The ideas of the proof are similar to Theorem [j] only that due to 
the added magnetic diffusion we need less information on the data. The main result 
of this section is that this decay is optimal. Specifically it is shown that the decay 
can not be uniform. 
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The proof we give below is formal. To make it rigorous, it can be applied to 
smooth approximations and then pass to the limit. The approximations could be 
constructed by retarded mollification as was done for the Navier-Stokes equations 
in [CKN82, Ler34 |. This construction if modified for the MHD equations will give 
suitable approximations which can be used to make our arguments rigorous. This 
arguments are standard and as such will be omitted. To see the construction of 
these approximations in detail we refer the reader to [HX05b|. 

5.1. Non-uniform decay. 

Theorem 2. Let n = 2, 3. If (u, B) is a weak solution to the MHD equations with 
S > 0, and data (u ,B ) <E (L 2 (R n )) , then 



lim 

t — >oc 



Mt)|l2 + ||fl(i)|| 



Proof. Without loss of generality suppose 6=1. The proof is based on similar 
arguments given in ]ORS97| ] for solutions to the Navier-Stokes equations with a 
forcing term. 

Let = e~\£> . As before split u into low and high frequency parts 



IK*)II 2 = l|fi(*)l 

\W)\\%= B(t) 



2 - ll<HI 2 + IK 1 - V)"ll2> 



< 
< 



ipB 



(l-<p)B 



Low frequency decay. We need to use Corollary (0) and Plancherel's identity. 



\<pu(t)\\ < 



-\t\ 2 (t-s) 



tpu(s) 



(<p 2 *u- Vu,e 2A{t ~ T) u) 



dr 



+ 2 f (^ 2 *5-VB !e 2A(i A) 

J S 



(5.1) 



dr = h 



and 



tpB(t) 



< 



l«| 2 (*-s) 



<pB(s) 



dr 



& 2 



VB,e 2A(t - T) 



dr = Ji + j r 2 + J3 



(5.2) 



It is immediate that the first terms I\ and J\ in (5.1), ( |5.2| ) tend to zero respectively 
as t goes to infinity. Hence it will be only necessary to show that the two integrals 
on the right hand side of each of the above equations tends to zero when t goes 
to infinity. Since all integrals can be estimated in a similar fashion we will only 
analyze the integrals I2, and ^3 corresponding to the velocity. 

Since (p 1 is a rapidly decreasing function, by the Hasdorff- Young, Holder, and 
Sobolev inequalities we have the following. 

When n = 2 



(ip 2 * u ■ Vu, e 2A(t ~ T) u) < (u-Vip 2 * 



2A(t- 



1^1 



i:i|Vu|| 2 



U, U) 

<C\\u 



< IMI: 



al|V«||2 



-2 + e 2A(t-r) Vu 
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and 



< 



< 



(B ■ V(p 2 * e 2A ^ T ^B, ; 



< \\B\ 



2 *e 2A ^~ T ^i 



|^ 2 |l ||B|i:i|Vu|| 2 <C||B|| a ||VB|| 2 ||V«|| 2 



<C(J|VB|| 2 + ||Vu|| 2 , 
For n = 3 

(<p 2 *u- Vu, e 2A ^u) | < ||^ 2 * u ■ Vu|| a ||«|| 2 < C ||^ 2 || 6/5 ||« • Vu|| 3/2 |M| 2 
<C|| U || 6 ||V U || 2 <C|| U || 2 ||V U || 2 . 

In the same fashion 



{Cp 2 *B- VB,e 2A{t - T) u) < ||<p 2 *B-VB|| 2 ||u|| 2 < C \\f> 2 \\ e/5 \\B ■ VB\\ 

<C\\B\\ 6 \\Vu\\ 2 <C\\u\\ 2 \\VB\\ 2 2 . 
Hence integrating over (s,t) yields 



3/2 H U ll2 



h+h<C j ||VB|| 2 + ||V U || 2 dr. 



Thus 

lim I 2 + h< hm lim C [ || VB|| 2 + || V?i|| 2 dr = 0. 

t — >oo s — >oo t — >oc J 

In the same manner it follows that the lim^oo J 2 + J3 = 0. Hence 



ipB 



0. 



(5.3) 



lim ||yu|| 2 = 0, lim 

t — >oo t — >oo 

High frequency decay. To estimate the high frequency part we will use Fourier 
Splitting [3ch85|. We now use Corollary (^|) and an equivalent version of this 
corollary for the magnetic field. Choose ip = 1 — = 1 — ip (note that in this 

case 4> is independent of time). Let x(t) = {£ <E K™ : |£| < G(t)}, then 



E(t) 



||(1-^)|| 2 + (l-cp)B(t) 



+ / E'{t) / |(1 -v)u(t)\ 2 d£dT+ E'{t) / (l-^B(r) 



<E{s) ||(l-^( s )|| 2 + (l-cp)B(s) 



d£dr 



x(t) 



x(t) 



E'(r) 



R*\x(t) 



|(l-^(r)| 2 d£dr+ / E'(t) 



m\ x (t) 



-2 E(t)U{1-^)u{t)\\1 dr-2 E(r) £(1-<p)B(t) 



(l-<p)B{r) 



dr 



d£dr 



E{t)(u ■ Vu, (1 - Lp) 2 u(r))dT + 2 / E(t)(B ■ VB, (1 - Lp) 2 u{r))dT 



2 E{t)(B ■ Vu, {l-ip) 2 B(T))dT + 2 / E(t)(u- VB,(1- tp) 2 B{T))dT. 



Suppose 



E(t) = (1 +t) a , and G 2 (t) = 



2(l+t) ' 



(5.4) 
(5.5) 
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with a > 3. With this choice E'(t) - 2E(t)G 2 (t) = 0, proceeding as in Theorem 
will yield 



E'(r) 



R?\ X (t) 



|(1 ~p)u{r)\ 2 d^dr-2 I E(t)U(1-<p)u(t)\\1 dr 



+ E'(t) (1-9)B(t) <%dT-2 E{t) Z(1-<p)B(t) 



R*\x(t) 



By (|5.6|), Equation (|5.4[) can be reduced to 



dr < 0. 
(5.6) 



E(t) 



\\{l-y)u{t)\\l+ (l-cp)B(t) 



<E(s) 11(1-^00^+ (\-*)B{8) 



f E'(t) f \{1- V )u{r)\ 2 d^dr+ [ E'(r) f 

J X (t) Js J X (t) 



t ft 



(5.7) 

d£dr 
(5.8) 



2/ E{t)(u-Vu, {1- ipyu{T))dT + 2 J E{r){B-^B,{l-ipYu{r))dT 

(5.9) 

t ft 



2 / E{t){B -\7u,{l- LpYB{T))dT + 2 I E(t){u • VS, (1 - ip) 2 B(T))dT. 

(5.10) 



We now will bound the terms in (5.8). Observing that |1 — tp\ < |£| if |£| < 1, we 
have 



l(i-^)l 2 H 2 



X(r) 



B 



d£ < CG(r)* I 

X(r) 



B 



<c(||« |l3 + l|Bo||2)(l + r)- 



(5.11) 



We now analyze (5.9), (5.1C) together. For this note first that (1 — ip) 2 = 1 + 6*, 
where 6 = —2tp + ip 2 , hence by the definition of ip, the function 9 is a rapidly 
decreasing function. 

Since (u ■ Vit, u) = (u- WB, B) = 0, and (B^VB, it) - (iTwB, B) = 0, it follows 



that the four last terms of the right hand side of (5.4) can be expressed as 



-2 / E(t)(u- Vu,6u{T))dT + 2 / E(t)(B ■ VB,6u{ T ))dT 

J S J S 

-2 j E(T)(B^Vu,6B(T))dT + 2 J E(t)(u^~VB, 6B(r)}dr = K± + K2 + -K3 + K4. 

The estimates of K^s are all very similar. Hence we only estimate K\ and state 
estimates are for Ki for i = 2, 3, 4. 
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For n = 2 



Kt= E{t)\u-Vu,6 *u(r)\dT = / E(t) I {u ■ * Vit(r), u(t)) I dr 



<C(\\uof 2 + \\B \\t 



EWpWJuWlWVuWzdT <C / E(T)\\u\\ 2 \\Vur 2 dr 

J S 

E(r)\\Wu\\ 2 2 dr, 
Same type of computations yields that 

Ki < C J E(t) (||Vu||2 + HVBH2) dr, i = 2,3,4. 
For n = 3 we also only estimate K\ 

E(T)\(u-Vu,6*u(T))\dT < I E(t) \\0\\ 6/5 \\u- Vu|| 3 \\u\\ 2 dr 



<C I E{T)\\u\\ & \\Vu\\ 2 dT <C I E(t) \\Vu\\ldr. 

J s J s 



The same estimates yield 



Ki<C J E(t) (\\Vu\\l + IIVBH2) dr, i = 2,3,4. 

Combining the estimates (5 1 7)-(^]l0), and the estimates for the K^s yields after 
division by E(t) 



\\(l-^)u(t)\\ 2 2 + (l-<p)B{t) 



E(s) 



2 - E(t) 



\\(1-V)u(s)\\ 



E(t) 
C 

W) 



E'(r) / \l-^[\uf + 

Jx(r) 



D 



S(r)(||Vu(r)||^ + ||VB(r)||^)dr. 



Since 1 — tp < 1, combining the last equation with (5.11), recalling the definition of 
E(t) and G{t) in (|5.5|), and since a > 3 we have 

2^ 



lim ||(l-¥>)«(t)lla+ (I-VW) 



< lim 



1 + s 



t^oo V 1 + t 

t 



ll^o || 2 + ll-Bolla 



+ 0(11^ + 11^)^(^^/(1+*) 
+ t lin i^)/ ts(T) (l" V ^ + l' V ^) rfT 



dr 



\\Vu\\l + \\VB\\l ) dr 



l*I| 2 -r || v ^n 2 ^ 

Letting s — > 00 on the right hand side yields 



lim [\\(l-ip)u(t)\\t+ (l-<p)B(t) 



= 
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Combining (5^3) with the last limit gives 



lim 

t — >oo 



\B\\ 2 =0. 



As stated in the beginning to make this proof rigorous, the formal proof has to be 
applied to the approximating solutions described at the beginning of the section, 
and then pass to the limit. This procedure is standard and as such is omitted. This 
completes the proof of the theorem. □ 

5.2. Lack of uniformity. In this subsection it is shown that solutions with data 
(uo, Bo) in (L 2 (K™)) cannot be expected to decay at a uniform rate, in the sense 
that for each sphere in L 2 (R") of radius (3, there is a point on the sphere so that 
the solution with such data will decay arbitrarily slow. In other words, given a time 
T > 0, and s > 0, there exists data uo with ||wo!| 2 = /? so that 

\HT)\\ 2 



> 1 -£. 



(5.12) 



Similarly for B given a time T > 0, and e > 0, there exists data Bq with H-B0II2 = 7 
such that 

\\B(T)\\ 2 



\B t 



> 1 



M 2 



(5.13) 



Theorem 3. There exist no functions G(t,(3) and if (£,7) with the following two 
properties. If (u, B) is a solution to equations ( jj.l[ )-( [77^ with S > and data 
(w , -Bo) € (L 2 {R n )) 2 , n = 2,3, then 

i) \\u(t)\\ 2 < G(t, \\u \\ 2 ), and \\B(t)\\ 2 < H(t, \\B \\ 2 ), 

ii) lim t ^o Git, [3) = 0, for all (3 > 0, and lim t _» H(t, 7) = 0, for all 7 > 0. 



Proof. As pointed in [3ch86| this lack of uniformity is already present at the level 
of the heat equation. 

The proof of Proposition 2.1 in [Sch8£] has a gap that will be taken care of in 
our present work. 

Notice that it suffices to show that ( 5.12| ) and (5.13) hold. The plan is to choose 
as initial data a family {uq , By } which satisfy 

and ||B a || 2 = j|i? || 2 . (5.14) 



-'o 1 1 2 



"0 



Write the solutions in integral form 
u a {x,t) =K*u%- 



[ K(x -y,t-s)* (u a Vu a - B a VB a + Vp Q )ds, 
Jo 

B a (x, t) = K*B%- [ K(x-y,t-s)* (u a VB a - B a Vu a )ds, 
Jo 



then 

|K(x,t)|| 2 > \\K*u%\\ 2 



[ \\K{x -y,t-s)* (u a Vu a - B a VB a + Vp Q ) 
Jo 

\\B a (x,t)\\2 > \\K * B0I2 - / \\K(x -y,t-s)* (u a VB a - B a Vu a )\\ 2 ds. 

Jo 



ds, 



(5.15) 
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The choice of (uq , Bq) will be the following 

vft(x) = a n / 2 u (x), Bg(x) = a n / 2 B (x), for n = 2,3. 

Then it is easy to see that Uq , Bq are invariant under a in L 2 (i.e., Equation ( |5.14[ ) 
is satisfied). Hence it is only necessary to show 



I K * u 



o 



"'0 



> 1 — e, and 



\K * B, 



o 



I -Bo 



>. 



and 



/ ||X(.t - y,t - s) * (u a Vu a - B a VB a + Vp a )\\ 2 ds < e, 
Jo 

[ \\K(x-y,t-s) * (u a VB a - B a Vu a )\\ 2 ds < e, 
Jo 

for e sufficiently small. 

We also note that the data (uq,Bq) will yield for the linear part a self-similar 
solution, that is 



K * u%{x, t) = a n/2 u{ax, a 2 t), 
K*B%(x,t) = a n/2 B(ax,a 2 t), 



(5.16) 



Equalities ( 5.1 6| ) follow by uniqueness of the heat equation and since a n / 2 u(ax, a 2 t) 
and K*v,q(x, t) are solutions to the heat equation Uq. In the same way cW 2 B(ax, a 2 t) 
and K * Bq (x,t) are solutions of the heat equation with data Bq. We only will 
show the proof for the energy of the velocity since the proof for the energy of the 
magnetic field is almost identical. 

We show first that Equation ( p. 12 ) holds for the data Uq. 



\u a \ 2 dx = a 1 



\u{ax,a 2 i)\ dx 



\u(y,a 2 t)\ dy 



R" 



d£ 



-2|?| 2 Q 2 t 



Next, by Lebesgue Dominated Convergence theorem it follows that for each fixed 
t, 



lim 



-2\£\ 2 a 2 t 



fio(OI 



|«o| d£. 



Hence 



lim B# = 1, 
||«o(-)lla 



Now it is necessary to show 
lim 



im f \\K(x-y,t- s)* (u a Vu a ~ B a VB a +Vp a )\\ 2 ds = 0, 
Jo 



(5.17) 



(5.18) 
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Hence we analyze 

rT 



I \\K(x-y,t-s) * (u a Vu a - B a VB a + V P a )\\ 2 ds 
Jo 



< [ \\VK(t-s) 
Jo 



\ 2 (\\u a u a \\ 2 + \\B a B"\L + \n)ds 



\ T - \ W J 3 || 3 + W x B a \\ 2 |, 



(5.19) 

Here we used that 

\\p a \\ 2 <c(\\u a u a \\ 2 + \\B a B a \C 

We can suppose then that we have chosen Uq and Bq to be in H 1 . Now observe 
that for n = 3 

\\u a u a \\ 2 < C*|jVu Q |j 2 /2 , \\B a B a \\ < C||VB"||2 /2 

and for n = 2 

\\u a u a \\ 2 < C||V U Q || 2 , \\B a B a \\ 2 <C\\VB a \\ 2 . 

Also 

||V«|| a + ||VS|| a < C{\\Vu\\ 2 2 + \\VB\\ 2 2 ) 1/2 and 

\\Vu\\l /2 + \\VB\\ 3 2 /2 <C(\\Vu\\l + \\VBf 2 f 4 . 

Thus to bound the right hand side of ( 5.19 ) we need to estimate 

tp(t) = \\Vu a f 2 + \\VB<*\\ 2 2 . 

Note first that if we choose (u ,B ) G Hi x Hi will yield that ||Vu Q °|| 2 , ||VB Q °|j 2 
are arbitrarily small if a <C 1 . This follows since 

<p(0)=<po = J \Vu«\ 2 dx + J \VB%\ 2 dx 

a A \\7u {ax)\ 2 dx + [ a 4 |V B (ax)\ 2 dx 



= a 2 \Vu Q (x)\ 2 dx + J \VB (x)\ 2 dx 
= Co?. 

In order to estimate ||Vm"|| 2 and ||Vi? Q || 2 we use Prodi's inequality We consider 
two cases: 

Case 1: (n = 2). By Prodi, 



*(f)<Cjfrf.)*, 



which implies 

p(t) < ^e c ^M ds < ^ e c (IKIl2+ll B oll2) < <fi 
< Ca 2 . 
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Hence, 



/ \\K(x-y,t-s)* (u a Vu a - B a VB a +Vp a )\\ 2 ds < CT 1/2 a 2 , 
Jo 



(5.20) 



and thus (5.1q) follows when n = 2. 



By 5.18 



and (5.17) the theorem follows for the velocity in two dimensions. 
Case 2: (n = 3). By Prodi, 



Integrating, 



11 r 

7"T < C / cp(s) ds 

<^o <P(t) J 



<C(\\u \\l + \\B \\l)=C. 



Solving for tp(t) we get 
¥>(*) < 



Wo 



< 2Ca 2 , 



1 - C <p 1 - Ca 2 
where we have chosen a very small, say Ca 2 < i, to make the last inequality 



true. Hence, for the case n = 3, the expression (5.18) is negligible too 



Combining (5.18) and (5.17) yields the conclusion of the theorem for the 
velocity in three dimensions. 

The estimate for the magnetic field in 2 and 3 dimensions follows in an 
analogous fashion. This completes the proof of the theorem. 

□ 

6. Kato's estimates 
In this section we show that by some simple modification Kato's pioneering work 



on L p decay for Navier-Stokes equations [|Kat8J holds for the MHD equations with 
magnetic diffusion. The main difference is that his approximating solutions will be 
replaced by the corresponding ones from MHD. Thus, rewrite the MHD equations 
as follows 

d t v + Av + F{v) = 0, 

where A = (A U A 2 ) = -P(A,<5A), v = (u,B), F(v) = F(v,v) = {F 1 (v,v),F 2 (v,v)), 
and 

Fi (u, w, B, D) = P(uWw) - P(BWD) 
F 2 (u,w,B,D) = P(uVB) - P(BVu). 

Here P is the orthogonal projection of L 2 onto the subspace PL 2 , which denotes 
the collection of divergence-free elements of L 2 . 
Following Kato 

u n +i = u (t) + Gi{u n ,B n ), 

Bn+l = Bo(t) + G2{u n , B n ), 
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where 

Gi(u,B) = - [ e-^-^F^uislBis^ds 
Jo 

and 

(u (;t),B (;t)) = (e-^uo^.e-^'Bo^)) . 

Using these expressions appropriately in theorems 1,2, and 3 in Kato's paper flKat84| 
will yield the same results for the MHD equations. 

We want to show how these results can be used to extend the decay results for 
the MHD equations in two dimensions when combined with our L 2 results. We 
recall, for easy reference, Kato's first two theorems. In this case u stands for the 
solution to the Navier-Stokes equations 

Ut + u ■ Vw + Vp = Ait, div u = 0. 

Theorem (Kato 1). Let the initial data a € PL m . Then there is T > and a 
unique solution u such that 

t (i- m / g )/2 ue BC([0,T);PL q ) form< q < oo, (1.1) 

t i-m/ 3 , Vit g BC([0, T);PL q ) for m < q < oo, (1.1') 

both with values zero at t = except for q = m in (1.1), in which u(0) = a. 
Moreover, u has the additional property 

u e L r ((0,T 1 );PL q ) withl/r=(l-m/q)/2, m < q < m 2 /(m - 2), (1.2) 

with some < T\ <T. 

Theorem (Kato 2). There is A > such that if \\a\\ m < A, then the solution u in 
Theorem (1) is global, i.e. we may take T = T\ = oo. In particular, ||it(£)|| 9 decays 
like t^ 1 "™/?)/ 2 as t — > oo, including q — oo, and ||Vu(i)|| decays like fr~ ( 1_m / 2 «) ; 
including q = m. 

As stated before usingg (u„+i, B n+ i) as defined above and following Kato's proof 
with straightforward modifications yields 

Theorem 4. Let n = 2, 3. Suppose (u , B ) 6 (PL P n PL n (R n )) 2 , where 1 < p < 
n. There exists Ai > such that if ||mo||„ < Ai and ||i?o||„ < Ai, then the solution 
to the MHD equations with 8 > is global and for any finite q > p 

t {n/p - n/q)/2 (u,B) and t {n/p - n/q+1)/2 (Vu,VB) G BC ([1, oo\; PL q f . (6.1) 

Combining the results of Theorems [l] and ^ in the two dimensional case yields 
the following improved decay for the solutions to the MHD equations with 6 > 0. 

Corollary 5. There is A > such that for ||uo|| 2 < A the global solution of the 
equation (6.1) for q > to, and for 2 < r < q 

Urn t^ ||(«,B)|L = 

t— >oo 

Proof. It follows by interpolating L r between L 2 and L q and using the decay rates 
of the solutions corresponding to those Sobolev spaces. □ 



NON-UNIFORM DECAY OF MHD EQUATIONS 



19 



7. Appendix 



Proposition 3. Let (u, B) be a solution to the MHD equations ( l.i. ). Assume 
the initial data mq, Bq is in L 1 (]R 3 ) n L 2 (R 3 ). Then 



\u{t)\ <C 1 + 



1 

w 

where C is a constant. 

Proof. We start by taking the Fourier transform of Equation (|l.2|) 



Let us define 



Ut + u- Vu — B ■ VB + Vp = - |£| u. 



H(£, t) = u-Vu-B-VB + Vp. 



Then u t + |£| u = ~H(£, t), and this equation can be integrated using the method 
of integrating factor to get 

u(t) = M(0)e-l«l 2 ' - / H(^,s)e-^ 2 ^ ds. 



o 



Then, 



\u(t)\ < \u \ 



[ |fr(£,*)|e-l*l a <*-"> ds. 
Jo 



To bound |if(£,s)|, let us first bound \Vp\. For this, let us take the divergence 
operator in Equation (1.2) which yields 

d 



It follows that 



ivpi = iei m < £ ^(i^ fe i + i^i) 



lei 



< ieiX)(l jBijB *i + i« , ' t **D 



Hence 



\H(Z >s )\<C\Z\Y;(\uiuK\ + \BiB*\) 

<C\t\{\\uo\\l + \\B \\l) <C\£\. 
Since |uo| < || u o|li = C, it follows that 

1*1 < l«ol + C |£| / e-l«'«-->ds <C + £(l- e-lfl' 1 



which finishes the proof. 



□ 
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